function [ gm ] = geomean( block )
    [n, k] = size(block);
    gm = rand(1,k);
    gm(1,1:1:k) = prod(block(:,1:1:k)).^(1/n);
end

